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We use the scalar model with quartic interaction to illustrate how a nonperturbative variational 
technique combined with renormalization group (RG) properties efBciently resums perturbative ex¬ 
pansions in thermal field theories. The resulting convergence and scale dependence of optimized 
thermodynamical quantities, here illustrated up to two-loop order, are drastically improved as com¬ 
pared to standard perturbative expansions, as well as to other related methods such as the screened 
perturbation or (resummed) hard-thermal-loop perturbation, that miss RG invariance as we explain. 

Being very general and easy to implement, our method is a potential analytical alternative to deal 
with the phase transitions of field theories such as thermal QCD. 
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At sufficiently high temperature or density, one could 
naively hope that the asymptotic freedom property of 
quantum chromodynamics (QCD) would give a reliable 
perturbation theory (PT) handle on the quark-gluon 
plasma physics. However, it is well-known that severe in¬ 
frared divergences unavoidably spoil a standard PT ap¬ 
proach in thermal QCD, and generically also for other 
thermal field theories, such that PT gives poorly conver¬ 
gent and furthermore badly scale-dependent results at 
successive orders (see e.g. [1] for reviews). Nowadays 
the development of powerful computers and numerical 
techniques offer the possibility to solve these nonpertur¬ 
bative problems in silico, employing lattice field theory 
(TFT). So far LFT has been very successful in the de¬ 
scription of the QCD phase transitions at finite temper¬ 
atures and near vanishing baryonic densities, with re¬ 
sults [2] which can be directly used for interpreting the 
experimental output from heavy ion experiments envis¬ 
aged to scan over this particular region of the phase dia¬ 
gram. However, the well-known numerical sign problem 
[3], which plagues this method when one considers the 
possibility of a particle-antiparticle asymmetry (signaled 
by a finite chemical potential), prevents LFT to be suc¬ 
cessfully used to describe compressed baryonic matter. 
Therefore, at the present stage, one cannot rely on LFT 
to describe the physics of compact stellar objects nor to 
explore the complete QCD phase diagram. In parallel 
over the last decades many efforts have been devoted to 
try to understand more analytically the bad convergence 
generically observed for thermal PT, even for moderate 
coupling values. Typically the dynamical generation of 
a thermal screening mass mu ^ \/\T influences the rel¬ 
evant expansion of thermodynamical quantities, such as 
the pressure, involving powers of vT rather than only A. 
Accordingly the predictions are, a priori, less convergent 
than for the T = 0 case. A plethora of nonperturbative 
approximations attempting to resum thermal perturba¬ 
tive expansions have been developed and refined over the 
years [1, 4-6]. The so-called optimized perturbation the¬ 
ory (OPT) is a variational approach in which a related 


solvable case is rewritten in terms of an unphysical pa¬ 
rameter, allowing for optimized nonperturbative results. 
In the past decades this strategy has been recycled, ap¬ 
pearing under different names [7-9]. At each successive 
order of such a modified perturbative expansion, the ar¬ 
bitrary variational mass is fixed by a stationary condi¬ 
tion. This strategy has already been used in a variety of 
different physical situations, including e.g. the determi¬ 
nation of the critical temperature for homogeneous Bose 
gases [10, II], the phase diagram of magnetized planar 
fermionic systems[I2, 13], and the evaluation of quark 
susceptibilities within effective QCD inspired models [14]. 
The development of a similar method, known as screened 
perturbation theory (SPT)[I5] or its version tailored to 
treat thermal gauge theories[16], hard-thermal-loop (re¬ 
summed) perturbation theory (HTLpt)[17], have been 
pushed to three-loop perturbative order [18-21]. Given 
the inherent technical difficulties of the (three loop) eval¬ 
uation of the QCD pressure for the case of hot and 
dense quark matter, the recent results in [21] repre¬ 
sent an impressive achievement. Moreover their agree¬ 
ment with LFT simulations is quite remarkable down to 
about twice the critical temperature, for the scale choice 
/i = 2ttT in the MS renormalization scheme. However, 
the SPT/HTLpt presents several shortcoming overshad¬ 
owing its potential as a reliable nonperturbative alterna¬ 
tive to LFT. Perhaps the most embarrassing issue is the 
strongly enhanced scale dependence displayed at increas¬ 
ing two- and three-loop orders, at odds with intuitive 
expectations: at three-loop order even moderate scale 
variations dramatically affect thermodynamical quanti¬ 
ties by relative 0(1) variations [18, 20, 21]. Another is¬ 
sue with the standard variational methods such as OPT, 
SPT or HTLpt is that beyond lowest orders, optimization 
gives more and more solutions, with unphysical complex¬ 
valued ones, often leading to use alternative prescriptions 
like replacing the variational mass with a purely pertur¬ 
bative mass [20], therefore loosing valuable nonperturba¬ 
tive information. 

Recently, the OPT method at vanishing temperatures 
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and densities has been consistently combined with renor¬ 
malization group (RG) properties[22-24]. The resulting 
RGOPT gives stable and precise results for the Gross- 
Neveu mass-gap [22], and new independent determina¬ 
tions [24] of the basic QGD scale (A^) and related cou¬ 
pling as, or the quark condensate [25]. Moreover, of¬ 
ten unique and real optimization solutions can be ob¬ 
tained [24], by matching those solutions to the RG behav¬ 
ior for small couplings; and by using appropriate renor¬ 
malization scheme changes. 

Here, we take an important step forward, by showing that 
the RGOPT is also compatible with the introduction of 
control parameters such as the temperature. To illus¬ 
trate how to implement the procedure we have chosen a 
simple, yet versatile model so that one can easily grasp 
the basic ideas and follow the main steps when perform¬ 
ing a particular application. More detailed results and 
formulas are given elsewhere [26]. Aside from purely cal¬ 
culation difficulties, the method described in this Letter 
can be directly extended to a large class of models. 

We thus start by considering the Lagrangian for one 
neutral scalar field with a quartic interaction, 

^ == , ( 1 ) 

where we have introduced a generic mass term m. The 
textbook result for the two loop free energy (equivalently 
minus the pressure) is [18, 27] 



ln(.;2+a;2) + _ 



+-^‘, ( 2 ) 


where in the imaginary time formalism = 2'KTn 
(n = 0,1, • • •) represents the bosonic Matsubara frequen¬ 
cies and LOp = p'^ + m? is the dispersion relation. The 
sum-integral in (2) as usual represents the sum over Mat¬ 
subara frequencies and remaining integration with mea¬ 
sure using dimensional regularization to 

perform the integral. The one-loop part of (2) is 


(4^)2 J-p = - 


m 


2 ' 
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where /x is the arbitrary renormalization scale in the MS 
renormalization scheme, and = m^/(4e) represents 
the vacuum energy counterterm[18]. We can already ad¬ 
dress a crucial point by considering the one-loop part free 
energy (3). It is a trivial matter to check that the renor¬ 
malized result spoils perturbative RG invariance. Acting 
on Eq. (3) with the standard RG operator: 


d d d d 


in (4), those being at least of next order (0(A). This is 
a manifestation of the fact that perturbative RG invari¬ 
ance generally occurs from cancellations between terms 
from the RG equation at order and the explicit /i 
dependence at the next order (our normalization is 

/3(A) = dA/dln/x = 5oA^-|-6iA^-|-for the/3 function and 

7 m(A) = dlnm/dln/i = qpA -I- yiA^ -!-••• for the anoma¬ 
lous mass dimension, with [28] (47r)^5o = 3; (47r)^7o = 
1/2; (47r)^5i = —17/3; (47r)^7i = —5/12). Neverthe¬ 
less, perturbative RG invariance can easily be restored 
by adding a finite vacuum energy term, £q, to the ac¬ 
tion without changing the dynamics. Although this term 
is usually ignored, minimally set to zero in the (ther¬ 
mal) literature [17, 18, 20], we stress that it is instru¬ 
mental for perturbative RG invariance to be achieved. 
Not surprisingly, we claim it largely explains the degrad¬ 
ing scale-dependence at higher orders in other similar 
resummation methods like SPT and HTLpt, which ig¬ 
nore those finite vacuum energy terms. The subtraction 
in MS-scheme is conveniently written as[24, 25, 29]: 

^o(A,m) = (5) 

fc >0 

where the coefficients Sk are perturbatively determined 
order by order from RG invariance. In the normaliza¬ 
tion of Eq. (3) we find sp = [2(5p — dqp)]”^ = 87r^, 
so that when augmented with £q the renormalized free 
energy from Eq. (3) is RG invariant at the one loop 
level. This can be carried out to higher orders, to give 
Si = —1, S 2 = (23-I-36C[3])/(4807r^), etc. Note that 
the apparently singular behavior for A —>■ 0 in (5) will 
actually disappear from the final optimized free energy. 
We stress that the previous construction, being only de¬ 
pendent on the renormalization procedure, does not de¬ 
pend on temperature-dependent parts: at arbitrary per¬ 
turbative orders the Sk coefficients can be determined 
from the T = 0 contributions only. This is indeed well- 
known, and the non-RG invariant remnant part defines 
the so-called vacuum energy anomalous dimension, that 
has been calculated even to five-loop order for the gen¬ 
eral 0{N) scalar model [30]. Our independent results for 
the Sk are fully consistent with [30]. A subtlety is that 
according to Eq. (5), Sk is strictly required for RG in¬ 
variance at order A^, but contributes at order A^“^. So 
at order A^ one may minimally choose to include only 
So,'' ‘ Sk, or more completely include s^+i ^ 0, thus in¬ 
corporating higher order RG information. 

One can now proceed to apply the RGOPT resumma¬ 
tion, by first performing on the RG-invariant free energy 
the substitution which appropriately modifies its pertur¬ 
bative expansion: 

—>■ (1 — (5)^“ ; A —>■ 5X, (6) 


and noting that the thermal contribution J^o{T) is scale 
independent, yields a remnant contribution: —m^/2, not 
compensated by lowest orders terms from /3(A) or jm(X) 


where now m is an arbitrary parameter, and the role of a 
is explained below. One then re-expands at successive or¬ 
ders 6^, setting d —>• 1 in the final results. This procedure 
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is consistent with renormalizability [31-33] and gauge in¬ 
variance [29], whenever the latter is relevant. The arbi¬ 
trary mass parameter m is then most conveniently fixed 
by a variational optimization prescription [9]: 

QTr(*) 

(m, A, d = = 0, (7) 

and m ^ 0 determines a nontrivial mass m(A) with non- 
perturbative A-dependence. 

In most previous OPT [7] (similarly SPT[15] and 
HTLpt[16]) applications, the linear d-expansion was 
used, i.e. assuming a = 1/2 in ( 6 ) mainly for simplicity 
and economy of parameters. However, to preserve RG in¬ 
variance after performing ( 6 ), a is uniquely fixed [24] by 
the universal (renormalization scheme independent) first 
order RG coefficients, as we recall below. Note, once 
combined with Eq. (7), the RG Eq. (4) takes a reduced 
massless form 

■ O "I 

M^+/3(A)^ -F^"^(m,A,d = l) = 0 , ( 8 ) 

so Eq. ( 8 ) with the OPT Eq. (7) completely set optimized 
m = rh and g = g “variational fixed-point” values. 
Gonsider the one-loop Eq. (3), at T = 0, augmented by 
£o = —{m'^/X)so, where as discussed above, sq = 87 r^. 
Performing ( 6 ), expanding to order consistently, and 
taking afterwards d —>■ 1 yields 

(47r)2 J-f = m* -^(1 - 4a) - (^^ + J In £ 2 ) ■ (9) 

Then to satisfy Eq. (8) implies a = 70/^0 = 1/6. At 
this one-loop order the RG Eq. (8) gives no further con¬ 
straints, but at higher orders it fixes an optimized cou¬ 
pling, and a = 70 /&o guarantees that among both RG 
and OPT solutions, at least one (often unique) is consis¬ 
tent with the T = 0 standard perturbative behavior [24] 
for A —>■ 0, i.e. infrared freedom in the present case: 
X{g <C m) ~ [60 \n(m/ ft)]~^. 

Switching on thermal effects, it is convenient to express 
our results in terms of the one-loop renormalized self¬ 
energy including all T-dependence, S/j, exphcitly[l, 18] 

Sfl= 7 oA m 2 (^ln^-l^+rVi(^) , ( 10 ) 

with the thermal integrals {t = p/T and x = m/T): 




4r[i/2] 

r[5/2-n] 




■ ( 11 ) 


Then noting that + Wp) = 2S/j/A, the ex¬ 

act solution of the one-loop OPT Eq. (7) takes the form 
of a self-consistent “gap” equation: 


= (47r)^ bo (12) 


which is exactly scale-invariant by construction. To il¬ 
lustrate this more explicitly, it is convenient to use the 
high-T expansion m/T = a; <C 1 of Jn{x), e.g. Jo{x) ~ 
167r‘*/45 —47r^x^/3-|-87rx^/3-|-x‘*(lnx/(47r) -I- 7 B — 3/4) -|- 
0{x^). This approximation is actually valid at the 0.1% 
level even for a; < 1 , sufficient for our purpose since the 
RGOPT one-loop solution fh/T always lies in this range. 
In this case the OPT Eq. (7) is a simple quadratic equa¬ 
tion for X, with the unique physical {x > 0) solution: 



with Lt = ln[/L( e'>'®/( 47 rT)]. We stress that the varia¬ 
tional mass (13) is unrelated to the physical screening 
mass [34]. The corresponding one-loop RGOPT pressure 
reads (Pq = (7r^/90)T‘* is the ideal gas pressure): 


, 15 „ 15 „ 45/1 

IT " Vm 



Eqs. (12)-(14) have clearly a nonperturbative dependence 
in A, and are exaetly scale-invariant, upon using for 
A = X{p) the “exact” (one-loop) running: 1/A(^') = 
1/A(^) — bolnp'/p,, then l/( 6 oA(^)) -I- Lt is explicitly 
/r-independent. Thus Eqs. (13) and (14) only depend 
on the single parameter 6 oA(/ro), where pLo is some ref¬ 
erence scale, typically po = SttT. This is a remarkable 
result, recalling that we started from (3) augmented by 
—m"‘(so/A) being RG invariant up to neglected higher 
order (!1(A), but not yet resummed, while (13), (14) are 
all-order RG invariant, showing the resummation effi¬ 
ciency after optimization. Eq. (14), perturbatively ex¬ 
panded, gives for the first few orders /Pq ~ 1—5a/4-|- 

5-\/6a^/^/3 -I- 5 (Lt — 6)a^/4 -I- where a = boX. 

Eqs. (12)-(14) reproduce exactly at arbitrary orders the 
0{N) scalar model large A-results {e.g. Eq. (5.7) of 
[35]), as can be checked upon identifying the correct 
large-TV bo = l/( 167 r^) value [35]. These results are 
also equivalent to those (at two-loop order) in [36], if 
replacing bo = 3 /( 167 r^) by bo/3, as argued in [36]. As 
we keep the correct bo, Eq. (14) differs from standard 
perturbative pressure by X{fj,o) Aperi(Lo)/3: this is 
not a problem, simply a different calibration as A is 
yet arbitrary since the model is not fully specified by 
any data fixing a physical input scale, /xq. Indeed, this 
apparent discrepancy disappears if expressing our re¬ 
sults in terms of the physical mass: to see it, we solve 
Eq. (7) now for A(m), replace it in (14), it gives simply: 

= l- 15 xV( 87 r 2 )-pl 5 xV( 87 r 3 )-pe)( 10 - 4 a; 6 ). But 
here x = m/T is arbitrary as we already used (7) to fix 
A(m). Now taking for m the physical screening mass [34] 

mf ~ (A/24)P^(1 — •\/ 6 A/( 47 r) H-), exactly reproduces 

the first two terms of the standard physical pressure [ 1 ]. 
Eq. (14) is plotted in Fig. 1, compared with stan¬ 
dard perturbative expansions at one- and two-loop orders 
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with their notoriously bad scale dependence [1]. Note 
that at one-loop, including non-minimally si 7^ 0 in 
(5) is actually equivalent to a simple scale redefinition, 
/r —^ = /re“^, in all our results above. 

The two-loop 0{6^) contribution to the free energy, for 
5 = 1, takes a compact form in terms of 'Sn in Eq. (10): 


^0 = 







(15) 


where £q^ = —m"‘/[l/(36oA) -I- si/3] from (5), and by 
abuse of notation the finite part of this already renor¬ 
malized expression is meant. The exact two-loop OPT 
and RG Eqs. (7) and (8) can be written compactly as 


/opt = (1 - 5 ^) + + S'p (5 - 3k) = 0, 


/rg = h 


ilk-s) A 


i/3(2)(A)52 = 0, (16) 


with h = (47r)“^, /3(^)(A) = 6oA^ -|-6iA^, and the reduced 
(dimensionless) self-energy S{m,fj,,T) = 'Z.h /We 
also have from Eq. (10): Ep = 9^2 (E_r) = A(5'-|-to^«S") = 
7 oA [ln(m^//x^) —J2 (w/'T)]- One may also solve the OPT 
and RG equations in the high-T expansion approxima¬ 
tion, which is excellent up to large (rescaled) coupling 
g = -^A/24 ^ 0{\) values and gives exactly solvable 
cubic and quartic algebraic equations respectively, with 
unique physical solutions {fh/T > 0 etc) easily identi¬ 
fiable. The resulting OPT and RG solutions for fh/T 
and P/Pq are consistent with Eqs. (13), (14) for the first 
two order terms perturbatively re-expanded, but contain 
appropriate modifications at higher orders (detailed ex¬ 
pressions are given elsewhere [26]). 

The exact two-loop pressure P/Pq obtained from the 
RG Eq. (4), as a function oi g = sJ\/2A, is plotted 
in Eig. 1, with scale dependence from exact two-loop 
running, compared with one-loop RGOPT and standard 
perturbative one- and two-loop pressure. The RGOPT 
improvement on convergence and scale dependence as 
compared to standard perturbative results is drastic, al¬ 
though a moderate residual scale dependence appears at 
two-loop, visible on the figure for (rescaled) coupling val¬ 
ues g > 0.6. This is not surprising since the construc¬ 
tion relies on a two-loop truncated basic free energy. At 
one-loop RGOPT the exact scale invariance is due to 
the peculiar form of the exact running coupling perfectly 
matching (13). At two-loop RGOPT the residual scale 
dependence reappears first at order A^: /S.P^QQpj,{g) ~ 
(0.075In^/^o — 1-92)(7®, i.e. one order higher than 
the normally expected A^ from standard RG properties. 
Moreover, including non-minimally S2 7^ 0 (thus catching 
a RG part of the three-loop contributions) modifies the 
perturbative pressure only at order A^, but slightly im¬ 
proves further the (nonperturbative) scale dependence, 
as intuitively expected and seen in Eig. 1. More remark¬ 
ably with S 2 ^ Q the two-loop pressure almost coincides 



FIG. 1: RGOPT P/Po{g = ■^A/24) at one- and two-loop ver¬ 
sus standard perturbative and two-loop SPT pressures with 
scale dependence ttT < g < 4nT. 


with the one-loop result up to relatively large g ^ 1. In 
Eig. 1 we also compare with the OPT/SPT two-loop 
result, i.e. discarding £0 in (5), taking a = 1/2 in (6), 
and using Eq. (7); and another prescription using instead 
the screening mass[34], similar to the QGD HTLpt pre¬ 
scription [20]. Note that the missing one-loop RG in¬ 
variance from unmatched m'^ln/x terms in (3) remains 
somewhat hidden at one- and two-loop thermal expan¬ 
sion order, since perturbatively m'^ ^ A^, explaining why 
it plainly resurfaces at three-loop A^ order in SPT[18] or 
similarly HTLpt[20]. In contrast the RGOPT scale de¬ 
pendence should further improve at higher orders: built 
on perturbative RG invariance at order k for arbitrary m, 
the mass gap will exhibit remnant scale dependence as 

^ Ar^(l-|-In g)), thus the dominant scale 

dependence in the free energy, coming from the leading 
term —som^/X, should be C>(A^+^). 

Finally we can combine the OPT and RG Eqs. (16) to ob¬ 
tain the full two-loop RGOPT solution, fixing fh/T and 
A = 24g^ for a given input scale g. We find for g = 2 ttT: 
fh/T ~ 0.912; g ~ 0.825; Prgopt/^o — 0.907, and the 
scale variation for ttT < g < A-kT is consistent with the 
one above shown. 

In conclusion, we have shown how resummations of ther¬ 
mal perturbative expansions based on a variational mass 
should be appropriately modified to restore perturbative 
RG invariance, missed by previous OPT/SPT,HTLpt 
analogous methods. The resulting RGOPT has a dif¬ 
ferent interpolation prescription, Eq. (6), uniquely dic¬ 
tated by universal first order RG coefficients a = 70/5o- 
The RG equation gives an alternative constraint to de¬ 
termine the nonperturbative variational mass and cou¬ 
pling, instead of using solely the optimization (7). The 
RGOPT pressure has exact one-loop RG/scale invari¬ 
ance, and a scale dependence and stability at two-loop 
order drastically reduced up to relatively large coupling 
values as compared with most other resummation ap¬ 
proaches. For thermal QCD we anticipate a similarly 
improved scale-dependence and stability from appropri¬ 
ate RGOPT adaptations of HTLpt. 
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